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Abstract 


In  this  report,  certain  properties  of  the  multicomponent  plasmas  present  in 
the  group-V  semimetals  As,  Sb,  and  Bi  have  been  derived.  Notable  among 
these  properties  is  anisotropy  of  the  effective  masses  of  both  electrons  and 
holes,  which  in  turn  leads  to  anisotropy  in  the  plasma  frequencies  of  these 
materials.  Because  the  systems  of  interest  are  particles  in  a  host  matrix 
rather  than  bulk  materials,  it  is  likely  that  this  anisotropy  in  the  plasma 
response  will  be  averaged  out  in  some  way;  however,  this  problem  must  be 
examined  in  detail  before  such  a  conclusion  is  warranted. 
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1.  Introduction 


The  electrodynamic  properties  of  the  semimetallic  group-V  elements  (ar¬ 
senic,  antimony,  and  bismuth)  are  complex  and  ill-understood  in  many 
ways.  This  is  primarily  due  to  the  complicated  structures  of  their  conduc¬ 
tion  and  valence  bands,  which  have  many  maxima  and  minima  for  the 
Fermi  surface  to  intersect,  resulting  in  multiple  "pockets"  of  electrons  and 
holes  populated  by  carriers  at  all  temperatures  [1].  According  to  the  Som- 
merfeld  theory  of  metals,  each  minimum  or  maximum  should  give  rise  to  a 
different  species  of  electron  or  hole,  so  that  the  electrodynamic  response  of 
the  material  resembles  that  of  a  multicomponent  plasma.  To  make  matters 
worse,  all  these  materials  crystallize  orthorhombically,  so  that  their  charge 
carriers  (both  electrons  and  holes)  have  anisotropic  effective  masses. 

When  nanometer-size  samples  of  these  materials  (so-called  quantum  dots) 
are  prepared,  the  resulting  systems  are  expected  to  behave  in  unusual  ways 
when  probed  by  electromagnetic  radiation.  A  particularly  easy  way  to  cre¬ 
ate  such  samples  was  discovered  some  10  years  ago,  when  researchers 
found  that  imder  certain  conditions  low-temperature  growth  of  GaAs  by 
molecular-beam  epitaxy  (MBE)  resulted  in  material  with  large  numbers  of 
ultramicroscopic  inclusions  of  As,  up  to  1  percent  of  the  host  volume  [2], 
which  arose  from  the  aggregation  of  GaAs  antisite  defects  (i.e.,  lattice  sites  at 
which  a  gallium  atom  was  replaced  by  an  As  atom).  Although  theories  re¬ 
garding  their  nature  are  still  problematic,  there  are  strong  indications  that 
these  inclusions,  which  are  roughly  spherical  in  shape,  consist  of  metal¬ 
lic  arsenic.  Similar  behavior  has  been  observed  when  InAs  and  GaSb  are 
grown  in  this  way. 

There  is  ample  reason  to  believe  that  these  materials  could  be  useful  in 
device  engineering.  However,  before  this  usefulness  can  be  explored,  the 
individual  nanoparticles  must  be  properly  modeled  with  regard  to  interac¬ 
tions  with  electromagnetic  waves.  In  this  report,  I  derive  an  appropriate  set 
of  constitutive  equations  that  can  be  used  with  the  Maxwell  equations  to 
implement  this  modeling. 
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2.  Boltzmann-Equation  Calculation  of  Fermi  Velocities  in  Group-V 
Semimetals 


The  most  effective  way  to  deal  with  metals  and  their  interactions  with 
electromagnetic  radiation  (light,  microwaves,  etc)  is  to  use  the  Boltzmann- 
Vlasov  equation  to  calculate  momentum  (A-space)  distributions  for  the 
metallic  charge  carriers  [3].  For  semimetals,  this  naturally  leads  to  a  gen¬ 
eral  multiband  set  of  Boltzmann  equations: 


where  t)  is  the  number  distribution  function  for  the  jth  carrier, 

Va  {k)  is  the  semiclassical  (hydrodynamic)  velocity  of  such  a  carrier  with 
momentum  k,  tj  is  its  particle-number-conserving  relaxation  time,  and 
0(a:,  t)  is  any  external  (electrostatic)  potential  applied  to  the  system.  The 
hydrodynamic  approximation  introduces  a  parameterization  of  the  solu¬ 
tions  to  these  equations  in  terms  of  macroscopic  variables  such  as  density, 
drift  velocity,  etc,  which  are  determined  self-consistently  as  moments  of  the 
equations.  If  the  particles  are  semiconducting  with  parabolic  bands,  we  can 
use  Maxwell-Boltzmann  functions  for  the  perturbed  carrier  distributions: 

/«'*£()  =  4A«)(£*)  e.p  <)}  .  (2) 

where 

G^^\k,  A)  =  y  (fc  -  k]^\x.  t)  -  ■  (3) 

Here  t)  is  the  hydrodynamic  density  of  carriers  in  the  jth  carrier 

pocket,  k'Y)  {x.  t)  is  the  hydrodynamic  drift  momentum  of  carriers  in  this 
pocket,  k^j  locates  the  center  of  the  pocket  in  A  -space,  n;  “  ^  is  the  effective- 
mass  tensor  of  the  jth  carrier  type  with  principal  values  nii  and  mi  (as¬ 
suming  the  pockets  are  spheroidal),  ks  is  Boltzmann's  constant,  T^^\x,  t) 
is  the  hydrodynamic  temperature  of  the  jth  carrier  type,  and  Ti  is  Planck's 
constant.  The  corresponding  unperturbed  distributions  are 


I 

I 
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where  n|/^  is  the  equilibrium  density  of  carriers  in  the  jth  carrier  pocket, 
and  T  is  the  usual  lattice  temperature.  For  metallic  particles,  the  correct 
distributions  at  low  temperatures  are  perturbed  Fermi  distributions: 

{k,  x.t)  =  9  (epix,  t)  -  epj  -  {k,  x,  t))  ,  (5) 


where  €f(x,  t)  is  a  position-dependent  Fermi  energy  (imref)  measured  from 
some  constant  reference  (e.g.,  vacuum),  and  epj  marks  the  minimum  en¬ 
ergy  of  the  jth  pocket.  The  corresponding  equilibrium  distribution  is 


where  ep  is  the  equilibrium  Fermi  energy.  We  can  relate  the  density  Uq  to 
cf  by  introducing  the  density-of-states  mass  /i  =  [mfini]  and  the  quan¬ 
tity  k^j  =  2iiAeFj/h?,  where  Aepj  =  ep  -  epj  (i.e.,  the  Fermi  wave  vector 
of  pocket  j),  and  by  integrating  the  distribution  over  fc-space  in  the  usual 
way.  This  gives 


n 


u)  _  tii 
°  3^2  ■ 


(7) 


We  can  likewise  define  the  deviation  from  equilibrium  in  the  usual  way: 

{k,  X,  t)  =  /(J)  (fc,  f ,  t)  -  {k)  .  (8) 


By  integrating  only  over  fc-space,  we  obtain  the  perturbed  and  unperturbed 
hydrodynamic  particle  densities  and  their  difference: 


rS^\x,t)  =  — ^  [  d^kf^\k,x,t)  , 
\J2i7r j  J 


5n^^\x,t)  =  n^^\x,t)  —  nQ\  (9) 


Note  that  these  also  enter  into  the  number-conserving  collisional  relaxation 
term  and  hence  must  be  determined  self-consistently.  These  equations  cou¬ 
ple  to  the  Poisson  equation  in  the  usual  way: 

V^(/)  =  -  ^  ej5n^\x,  t)  ,  (10) 

^  i 

where  the  charge  ej  =  +1  for  holes,  —1  for  electrons. 

Let  us  linearize  the  Boltzmarm  equations  around  the  unperturbed 
distributions: 

+v'i\k)a„6p'  +  (^/o'Vi))  =  ~  .  (11) 
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where  we  neglect  all  spatial  dependences  to  zero  order.  To  solve  these  equa¬ 
tions,  we  use  the  (wave)  ansatz 


■'  j 

Then  the  equations  for  the  perturbed  system  become 


(12) 


-ioj  +  Sf^^\k)  +  iqaV^^\k)Sfjj^\k)  -  -^(/V?o0g)  ^ 


'jrfo’ik)  ,  (13) 


n. 


(i) 


which  are  easily  solved: 

sfi^'X)  = - ^ - -■-■■  ukk - ^ 

<^+  j-  -  Qa'i’a  {k) 

We  can  now  substitute  this  expression  back  into  the  density  expression 


= 


i2^r 


.^Hk)  d^k 


(15) 


to  get 


which  can  be  written 


d^k  + 


i  Sn 


0) 


_  f _ fl 

)3  7  ^  ,  ±  _ 


Tj  (27r)'^  J  q^vl^^k) 


d^k  ,  (16) 


(17) 


j  n,. 


The  first  term  in  this  expression  contains  the  vector  coefficient 

2 


C>'7  _ 


3,  !ife4^'(*-) 


(27r)'^  J  i  _  q^vli\k) 

while  the  second  term  contains  the  scalar  coefficient 

2 


(18) 


^2j 


fdh 

(2^)"  7  E'+7--go4'>(0 


(19) 


4 


Since  each  /c-space  pocket  in  the  semimetal  is  ellipsoidal,  the  hydrodynamic 
velocities  can  be  defined  in  terms  of  the  effective  mass  tensors: 


v^J\k)  =  n 


m 


(20) 


where  is  the  deviation  in  k  from  the  jth  band  minimum  (or  maximum). 
Then  the  scalar  term  becomes 


= 


(27r) 

and  the  vector  term  becomes 


Aj)( 
Jo  V 

f) 

^  +  :r  “ 

‘j 

^7 

a7 

CV0  _ 

- 


(27r)' 


/ 


0J  +  ^  -qafi  i 


a'j 


Using  the  Fermi  distributions 


m. 


-1 


gives 


('r  -  -  Y«« 


J  ap 


J  ajS 


(21) 


(22) 


(23) 


(24) 


Then  becomes 


Qt  -  -n 


(27r)^ 


^  K  \ 

/37  ^ 

CO  +  ^  —  qah 
‘3 

1 

1 _ 1 

where 


[mj  . 

We  can  now  solve  the  density  equation 


S{eF-e^j)  (25) 


(26) 


1  - 


—^2i. 

Pj 


Sn‘p'>  =  ^(g/30q) 


hj 


and  insert  the  results  into  Poisson's  equation  to  get 


1  - 


i  Q2j 


‘3  n, 


O') 

0  J 


-1 


c>/5 


(27) 


(28) 


Note  that  for  solid-state  systems.  Equation  (28)  can  only  be  true  if 


1  = - v 


1  - 


i  ^2j 

T,-  , 

J  Uq  J 


-1 


1 

h 


cv,/3 


(29) 
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where 


cx.a  f 

q;3^ii  =  -Ti 


{2tt) 


,/A; 


H3 

rrij  ^ 

a 

J 

Ud+ij- 

hqrj 

57^^ 

(30) 


The  roots  of  this  equation  define  various  dispersion  relations  for  plasma 
waves  ^^*^((7). 

Because  of  the  ^-function,  it  is  possible  to  evaluate  expression  (30)  analyti¬ 
cally.  Write 


m , 


a  3 


—  SaX^Xi/S^ifi  . 


(31) 


where  the  matrix  5  is  a  rotation  chosen  to  diagonalize  m-  ^ .  Then 


rx./?  -f- 


(27r 


fi(};3SgxDxi,S,y^:^- 


"t”  X, 


6  Affj  -  — ^  (32) 


where  Aep-j  =  ep  -  (pj.  Let  x^,  =  Qy  =  S„y  qy.  Then  the  Jacobian  for 

this  transformation  is  1,  and  so 


„  0,5  _  2  A,,3„  TlQxDxuXy  s:  (  A  .  TA 

q0'^\j  .3  j  (I  X  .  t  r\  o  ^  X^\\/x'iy 

(2py  J  u  + ^ -hQxDxyXy  \  2 


TiQxDxyXy  „  (  .  ry 


(33) 


Now,  since  D  is  diagonal,  we  can  easily  define  D  —  as  a  scaling 

transformation.  If  Xy,  Ay  =  Ti  Qy,  then 


m^rrti 


li/2  2  y  3  A„y„ 

^jdy - 


n- 


{2pYJ  -a;+^-A,t/,, 


(5  Aepj  -  —y^  ,  (34) 


where  the  factor  in  front  is  the  Jacobian  of  the  scaling  transformation.  In 
spherical  coordinates,  this  expression  becomes 


0  5  r  2 

^0%  =  [^tm 


ll/2 


^  dy  5  (Aepj  -  y  j  ^  d{cos  0) 


Ay  cos  9 


ix  +  Y  —  Ay  cos  9 


(35) 


where 


r  /  ^2  2  \  1 

1/2 

A  = 

z= 

rYqa  (mj^)  a  3 

\mt  mi  J 

(36) 


Here,  q±_  and  q.  are  components  of  q  perpendicular  and  parallel  to  the  c-axis 
of  the  particle  material  in  the  principal-axis  system  of  cylindrical  coordi¬ 
nates.  I  have  used  the  fact  that  as  a  scalar  |  A|  must  be  rotationally  invariant. 
Introducing  the  density-of-states  mass  y  =  [ni^n?/]  and  the  Fermi  wave 
vector  Jipj  —  2yA€  fj/h^  for  pocket  /  as  we  did  above,  we  obtain 


(7vj3  _ 

Q3^lj  — 


(37) 
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where 


2AeFj  (qa  {mj 


contains  the  tensor  nature  of  the  effective  mass  both  explicitly  and  implic¬ 
itly  through  the  Fermi  energy.  The  integral  is  trivial: 

{2  -  rin  \  =  (T)  .  (39) 


27r^7i 


r-1. 


The  quantity  92j  can  also  be  evaluated  explicitly.  Analysis  similar  to  that 
given  above  leads  to  the  double  integral 


^  1  r  2  11/2  /■°°  2 

^  I  y 


f  y  dye  I  -  -y  1  J_^  d(cos0)-— .  (40) 


Introducing  the  notation  yp  =  e  =  cos  0  as  before,  and  using 

the  properties  of  the  Fermi  function,  we  can  rewrite  the  above  as 


^22  = 


1  r  2  11/2  2  ,  1 

m^m/  /  ydy  d^ — —■ - 

>-  J  Jo  J-l  iO  +  ^  — 


Some  algebra  gives  the  following  expression  for  this  integral: 


(-+4) 


-/2(r) ,  (42) 


where  F  is  the  quantity  introduced  above. 

Let  us  solve  the  plasma  dispersion  relation  to  lowest  order  in  q,  i.e.,  the 
long-wavelength  limit.  To  do  so,  we  take  the  small-g  limit  of  our  explicit 
expressions  in  equations  (30)  and  (42),  which  corresponds  to  F  — >  00.  Some 
algebra  shows  that  the  function  /2(F)  defined  in  equation  (42)  goes  to  1  in 
this  limit,  so  to  lowest  order 


92j  = 


J  a;  + 


^  + 


The  same  limit  for  qy'^^j  gives 


„,cv/3  r  1 

~  I  3F2 


37r2)j  I  [u}-\-i/Tj) 


7 


But 


2fi,AeFjkfj 


^  4.  51 

nil  nil 


Stt-Y? 


=  -Ti 


2  kf 


5i  +  Ml 
L  ) 

37r^  I  (w  +  i/TjY 


=  -r>n!/U 


so  that  equation  (29)  becomes 


(45) 


-/  1 
+  :r 


-1 


JJ) 


2  Qo 


rn 


-1 


aii 


(46) 


Assuming  that  the  relaxation  times  are  the  same  for  all  the  pockets,  we 
obtain 

.2 


1-4 


1 


MM 


(.+i) 


Qa  {  Y. 


(i) 

0 


77} 


-1 


<13 


q3 


or,  introducing  the  unit  vector  q  =  q/ 

0 

1-^. 


UJ 


0. 


(47) 


(48) 


This  expression,  which  is  clearly  in  the  form  of  a  requirement  that  a 
"Drude"-type  conductivity  for  the  system  should  vanish,  defines  the 
"plasma  frequency"  of  the  multicomponent  system  as  the  sum  in  the  nu¬ 
merator.  In  general  this  sum  will  depend  on  the  direction  of  q. 

To  derive  an  analogue  of  the  linearized  Thomas-Fermi  equation  and  a  cor¬ 
responding  screening  length,  we  need  to  take  these  expansions  to  higher 
order  in  q.  Saving  next-order  terms  in  /]  (T)  and  /2(r)  gives 

/i(r)  «  ^  .  Mr) » 1  +  i  .  (49) 


Let  us  rewrite  equation  (29)  as 


1  = 


f<7‘ 


E 


IMl 

Tj  OJ  +  ^ 


p.k 


Fj 


Ti'^h 


h  (F) 


and  note  that 


/jFfj  _  3 

2  Aefj 


(50) 


(51) 
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Then  equation  (50)  becomes 


1  - 


ea: 


±MX 

to 


-1 


/i(r) , 


(52) 


where  is  the  Thomas-Fermi  screening  length  for  carriers  in  the  j  th  pocket; 


2  _  o  riQ  e 
^  2e  Ae pj 


(53) 


With  these  results,  we  can  easily  derive  the  expression 


i  =  E 

j 


1  + 


6  Acfi 


(-n)' 


2  i  ' 

1  +  o  r 

3  LU'T' j  . 


Qa 


m 


-1 


a/3 


(54) 


The  peculiar  fraction  (6/5)  that  occurs  in  this  expression  derives  from 
our  linearization  around  a  spatially  independent  zero-order  distribution. 
A  more  careful  linearization  is  needed  to  recover  the  true  hydrodynamic 
value  of  this  coefficient  and,  indeed,  to  recover  any  of  the  expressions  en¬ 
countered  in  classical  fluid  mechanics.  This  will  be  the  subject  of  a  subse¬ 
quent  report. 
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3.  Calculations  of  Sums  Over  /c-Space  Minima  for  High-Symmetry 
Band  Structures 


3.1  Simple  Cubic  Band  Structure 

A  simple  cubic  band  structure  is  appropriate  for  elemental  silicon  (see  fig. 
1).  In  order  to  calculate  the  effective  masses  for  this  method,  we  assume 
that  the  densities  in  the  carrier  pockets  are  the  same,  and  that  the  principal- 
axis  projections  are  xx,  f/y,  and  if.  Since  there  are  six  pockets,  each  holds 
1  /6  of  the  electrons. 

Let  us  write  the  effective  masses  as  follows: 


m*  ^’1  =  —xS:  +  —  {fjy  +  if)  . 

J  nil  m-t 

m'-  =  — yy  +  —  [xx  +  zz)  , 

J  mi  nif 

177^  ^  —  — zz  -] - (yy  +  xx)  . 

-I  nil 


Then  the  sum  is  simple: 


H  ^0  ^  ["'j  ^1  =  (—  +  — )  2  [.f.r  +  yij  +  if]  =  i/io  (—  +  — )  T  , 

^  L  J  J  6  \mi  miJ  '■  '3  \mi  nitJ 


where  1  is  the  identity  matrix.  This  shows  that  the  mass  required  is  just  the 
usual  isotropic  optical  mass. 


Figure  1.  Conduction 
band  minima  (pockets) 
of  silicon  in  fc-space. 
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3.2  (111)  Cubic  Pockets  (Tetrahedral  Band  Structure) 


A  tetrahedral  band  structure  is  appropriate  for  elemental  germanium  (see 
fig.  2).  Here  we  have  four  pockets,  at  the  vertices  of  a  tetrahedron.  The 
effective-mass  tensors  are 
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rrit 

[siSl 

+ 

iitij 

m(2)- 
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-2)' 

-1 

1  .  . 

- 6262 
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+ 
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mt 

(^S2S2 

-f- 

r^(3) 

-1 

m^~ 

-3)' 

1  .  . 
— 6363 

+ 

1 

(S3S3 

kh) 

L  J 

mi 

m-t 

/ 

-1 

m^~ 

-4)' 

-1 

1  .  . 
— 6464 

-t- 

1 

(S4S4 

+ 

*4^4) 

L 

- 

mi 

m 

\ 

/ 

and 


(57) 


where  the  vectors  ei,Si,ii  =  ei®  Si  mark  the  principal  axes  of  the  pocket 
ellipsoids: 


ei  = 

^{x  +  y  +  z)  , 

5i  = 

62  = 

j={-x  +  y-z)  , 

1  h  ^ 

(  h  — 

63  = 

j.{-x-y  +  z)  , 

)  53  = 

k  = 

64  = 

1 

1 

54  = 

j^{i  +  y)  , 

U  = 

^(-i-y  +  2i)  ; 

{—X  y  +  2z)  ; 
■^{x  +  y  +  21)  ; 
^  (x  -  y  +  2i)  . 


and 


(58) 


Figure  2.  Conduction 
band  minima  (pockets) 
of  germanium  in 
A:-space. 
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Then  the  dyadic  products  can  be  written  in  matrix  form  as  follows: 


1  1  1 
eiei  =  i  [ill 
1  1  1 


ese-i  =  4 


1 

-1 

1 

-1 

1 

-1 

1 

-1 

1 

1 

1 

-1 

1 

1 

-1 

-1 

-1 

1 

1 

-1 

-1 

-1 

1 

1 

-1 

1 

1 

•  54^4  = 


r  Uh  =  n 


hh  =  k 


•  hh  =  g 


^4^4  =  I 


(59) 


Performing  the  summation  once  more,  we  obtain 


m 


(1) 


3  m/ 
1 

mi 


/I 

1 

IN 

)+-■ 

fi  / 

1 

-1 

0 

1 

1 

1 

- 1 

-1 

1 

0 

Vi 

1 

1> 

/  'ffit 

in 

0 

0 

0 

(1 

1 

3 

h 

1 

i 

i 
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1 

1 

1 

2 

1 

3 

3 

3 

_| - 

3 

3 

3 

1 

1 

1  , 

nit 

1 

1 

2 

\3 

3 

3  / 

V  3 

3 

3  / 

m 


(2) 


-1 


1 

2 


m. 


(3) 


Smi 

J_ 

m,/ 


1 

0 


3 

3 

3 

3 

3 

3  ^ 

1 

1 

1 

1 

1 

2 

1 

3 

3 

3 

_| - 

3 

3 

3 

1 

1  Q 

1 

•3 

4 

771 1 

1 

l  o 

1 

2 

o  i 

(60) 
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lO  to 


and  so 


Each  pocket  holds  1/4  of  the  electrons.  Then  the  total  mass  sum  equals 


which  is  again  the  isotropic  optical  mass. 
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4.  Calculations  of  Sums  Over  /c-space  Minima  for  the  Band  Structures 
of  Group-V  Elements 


The  band  structures  of  the  group-V  elements  are  described  in  detail  by  Lin 
and  Falicov  [4]  and  Shapira  and  Williamson  [5].  For  the  conduction  band, 
the  pattern  of  ellipsoidal  pockets  shown  in  figure  3  is  obtained;  the  angle 
with  4;-axis  8c  is  material-dependent.  I  assume  here  that  the  valence  band 
has  a  similar  structure,  although  the  actual  Fermi  surface  of  As  is  a  single 
multiconnected  surface  rather  than  a  collection  of  pockets. 


Let 

a  =  cos  6c  ,  T  =  sin  6c  . 


(63) 


Then  the  hexagonal  symmetry  in  the  basal  plane  motivates  us  to  introduce 
the  quantities 

0  =  1/2,  d  -  V3/2  .  (64) 


This  allows  us  to  write  the  longitudinal-mass  principal  axis  vectors  as 
follows: 


ei  =  Tx  -t-  az  , 

€2  =  OiTx  +  f3ry  —  az  , 

63  =  —arx  +  Pry  +  az  , 

64  —  ~rx  —  az  , 

65  —  —arx  —  jSry  +  az  .  and 

66  =  arx  —  (irij  —  az  .  (65) 


Figure  3. 

Conduction-band  and 
valence-band  minima 
(pockets)  of  As  and  Sb 
in  A:-space.  Upper 
pockets  are  filled  with 
electrons,  lower  ones 
with  holes. 

(a)  Three-dimensional 
view,  and  (b)  view 
down  c-axis  (z-axis  in 
(a)). 
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Then  the  other  principal  axes  are 


eiei 

62^2  = 

6363 

6464  = 

hh  = 

^6^6  “ 


51  —  y  ,  ii  =  —ax  +  rz  ; 

52  =  —l^x  +  ay  .  h  =  aax  +  (3ay  +  tz  ; 

53  =  -(5x  —  ay  ,  is  =  —aax  +  Pay  —  rz 

54  =  —y  ,  i4  =  —ax  +  TZ  ; 

55  =  Px-  ay  ,  4  =  aax  +  Pay  +  tz  ; 

sq  =  Px  +  ay  ,  h  =  aax  —  Pay  +  tz  . 


and 


(66) 


As  before,  we  find  the  projection  operators 


,ar  0 

a3r‘^ 


0 

ar\ 

0 

0 

1 

0 

al3r^ 

—ara  \ 

-/3t<j  1 

—Pra 

/ 

—aiSr'^ 

—ara 

SlSl 


5252  = 


5353 


5454  — 


S555  = 


—ara 

—Pra 

V  0 

0 

a^T^ 

—apr"^ 

—ara  \ 

ap 

—afSr^ 

P^r^ 

Pra  1 

,  5656  =  OtP 

—ara 

Pra 

/ 

\  0 

0 

— a/3cr^  0^(7'^ 

arc  ~l3rcr 


The  dyadic  sums  are  then 


i=l 

6 

^  ^  ii  ii  ~ 

i=l 

6 

= 

i=l 


'  4a^T^+2T^  0 

0  4/32t2 

,  2(l-2a)T<T  0 

4a^a^+2a^ 

0 

— 2(1— 2q)to’ 

'4/32  0 

0  4a2+2  0 

0  0  0 


(67) 


,  and 


(68) 


Hence, 

6 


^  ^2 5^  iiii  — 


i=l 


'  4/3^+4a^cr^+2cr^  0  —2{l—2a)ra^ 

0  4a2-h4/?2a2+2  0 

—2(1— 2a)rcr  0  6r^ 


(69) 


Note  that  because  a  =  1  /  2,  this  tensor  is  diagonal!  Since  P^  -3/4,  the  above 
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becomes 


6  /  3{1+(t^)  0  0 

+  iiij  =  I  0  3(i+(t’)  0 

\  0  0  , 

Likewise,  we  have 

0  /3r2  0  0 

=  0  3r2  0 

i=y  \  0  0  6(7^ 

so  that  the  averaged  optical  effective-mass  tensor  is 


(70) 


(71) 


E 


71, 


U) 


ni' 


(  3  (1  -|-  <7 
0 
0 


0 

3(l-fcr^)  0 

0  6r“ 


0  \1 


0  0 

+  w  (1  +  0 

0  + 


/J 
\ 


(72) 


In  keeping  with  our  expectations,  the  effective  mass  tensor  is  somewhat 
anisotropic.  However,  there  is  an  angle  at  which  it  becomes  isotropic;  this 
is 

cos  200  =  —^ — .  (73) 

nit  — 

It  is  easy  to  show  that  the  reported  angles  for  the  group-V  materials  do  not 
satisfy  this  criterion,  so  that  anisotropy  will  be  a  real  issue  in  any  study  of 
their  plasma  response. 
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5.  Conclusions 


In  this  report,  certain  properties  of  the  multicomponent  plasmas  present  in 
the  group-V  semimetals  As,  Sb,  and  Bi  have  been  derived.  Notable  among 
these  properties  is  anisotropy  of  the  effective  masses  of  both  electrons  and 
holes,  which  in  turn  leads  to  anisotropy  in  the  plasma  frequencies  of  these 
materials.  Because  the  systems  of  interest  are  particles  in  a  host  matrix 
rather  than  bulk  materials,  it  is  likely  that  this  anisotropy  in  the  plasma 
response  will  be  averaged  out  in  some  way;  however,  this  problem  must  be 
examined  in  detail  before  such  a  conclusion  is  warranted. 
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